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Abstract 

We discusse a relativistic Hamiltonian for an n-body problem 
in which all the masses are equal and all spins take value 1/2. In 
the frame of reference in which the total momentum P = 0, the 
Foldy-Wouthuysen transformation is applies and the positive 
energy part of the Hamiltonian is separated. The Hamiltonian 
with unharmonic oscillator potential is applied to describe mass 
differences for charmonium and bottonium states. 
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1 Introduction 



In the development of non-relativistic quantum mechanics the contri- 
bution of Schrodinger was not only in writing a single particle wave 
equation that, for a Coulomb potential, leads to the spectrum of the 
hydrogen atom, but also that his formalism could be extended imme- 
diately to a system of n-particles in configuration space. Thus wave 
mechanics from the very beginning was able to deal with systems of 
many particles in interaction as happens for electron in atoms or nu- 
cleons in nuclei. 

The non-relativistic many body Hamiltonian, in the absence of 
spins, was composed of the sum of the kinetic energies Y^=\ip1l'^ m s) of 
the n particles together with an appropriate potential interaction be- 
tween them. Many refined mathematical methods [1] were developed 
to obtain the eigenvalues and eigenstates of these Hamiltonians. 

The corresponding development did not take place in the theory of 
many body problems in relativistic quantum mechanics. It is true that 
almost simultaneously to the development of non-relativistic quantum 
mechanics Dirac [2] introduced his famous relativistic wave equation 
for a single particle. The extension of this equation to many bodies 
found two main obstacles. One of them was the need to formulate 
the wave equation in a form invariant under the transformations of the 
Poincare group, thus guaranteing its relativistic character. In particular 
relativistic invariance presupposes a certain symmetry between space 
and time variables, which causes the appearance of multi-time variables 
in an n-body theory and thus is an aspect that we will correct without 
violating Poincare invariance. 

The other obstacle was that in relativity the relation between energy 
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E and momentum p of a particle of mass m is 



ztt2 2 2 i 2 4 

£/ = p c + to c 



(1) 



instead of the E = p 2 /2m of the non-relativistic case. 

Equation (1) implies that both positive and negative energies are 
possible as taking the square root we have 



This sign ambiguity makes itself felt immediately in the Dirac equa- 
tion and led him to propose that vacuum was actually a state in which 
all negative energy levels were filled by particles obeying Fermi statis- 
tics and thus unable to accept another particle of this type. From 
this point of view relativistic quantum mechanics turned to a field the- 
oretical description and a procedure similar to that of non-relativistic 
quantum mechanics for relativistic many body problems was essentially 
abandoned. 

There are many different formulations of the relativistic many-body 
problem [3]. Conventionally they can be subdivided to three classes. 
Namely, we can indicate manifestly covariant approaches based on the 
Bethe-Salpeter equation and its generalizations, theories with direct 
interaction and the mass shell constraints approach. 

The Bethe-Salpeter equation [4], [5] presents powerful and elegant 
tools for construction of various two body theories which are transpar- 
ently relativistic invariant. However, a generalization of this equation 
to n-body case seems to be too complicated if at all possible as far as 
practical applications are concerned., see, e.g., reference [6] for the case 
n — 3. 

A global receipt to overcome difficulties with multi time variables 
was proposed by Dirac [7] who proposed to use one time formulation for 
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n-body models but to ask for existence of realization of the Poincar al- 
gebra on the set of solutions of the equations of motion. The price paid 
for the absence of extra time variables was the absence of manifestly 
relativistic invariance. 

The mass shell constraints approach [8] shares with the Bethe- 
Salpeter equation manifestly relativistic invariance. This approach is 
much more easy to handle than ones based on the Bethe-Salpeter equa- 
tion. However, the problem of separation of extra time variables is by 
no means trivial and was effectively solved only for two- and three-body 
problems [9]. 

In this paper we want to return to the possibility of discussing rel- 
ativistic many body problems in a framework similar to that of the 
non-relativistic quantum mechanics. For this purpose we have to deal 
with two problems. 

a) How can we formulate a wave equation explicitly invariant un- 
der the Poincare group but which, in an appropriate frame of 
reference, involves only one time. 

b) Once objective (a) is achieved how can we separate the positive 
and negative parts of its solution as only the first one will be of 
physical interest. 

In Section 2 we deal with the first problem through a method we 
developed in previous publications and in section 3 with the second one 
employing a generalization of the Foldy-Wouthuysen (FW) transforma- 
tion. 
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2 A Poincare invariant n-body wave equa- 
tion which, in a particular frame of ref- 
erence, involves only one time 

We mentioned in the introduction how we can formulate the wave equa- 
tion for a n-body non-relativistic problem starting from the correspond- 
ing one body expression for the free particle. In the relativistic case we 
have also the wave equation of Dirac [2] for a single free particle given 
by 

d \ 
-i— + a-p + f3m\^ = (3) 

where we use the usual relativistic units h — c — 1, with p being the 
momentum three vector, m the mass of the particle and the matrices 
ex., (3 are defined as in reference [2]. 

If we have n particles of the same mass m, we add the index s = 
1, 2 ... n to all the variables and an obvious Poincare invariant n-body 
equation can be written in the form 

n 

E(7>^ + ^ = (4) 

s=0 

where repeated index \i are summed over the values \i = 0, 1, 2, 3 with 
Pos = —id/dt s , pj S = —id/dxj S ,j = 1,2,3 and the 7^ are matrices 
related to ck s and /3 S by [2] 

7° = Ps, ll = Ps^is, s = 1, . . . n, i = 1, 2, 3 (5) 

The 7^,p MS are respectively contravariant and covariant expressions 
[2] so that ^gPus is a Poincare scalar and thus Eq. (4) is certainly a 
Poincare invariant n particle wave equation, but it is not satisfactory 
because introduces n times through p 0s = —id/dt s , s — 1, 2, . . . , n. 

How can we find a formulation of many body problem, still invariant 
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under the Poincare group but, in an appropriate system of reference, 
involving only one time? 

We start by denoting by unit time like four vector which implies 
that there is a reference frame in which it takes the form 

(m m ) = (1,0,0,0) (6) 
With the help of the four vector (6) we can define the Lorentz scalars 

r=]I(7X), iWtX)-^ (7) 

r=l 

where (7^w M ) _1 eliminates the corresponding term in T so that T s is 
still in product form. 

Instead of Eq. (4) we propose now the following Lorentz invariant 
one [10] 

]rr s ( 7 £V + m)v> = o. (8) 

s=l 

We now introduce the total energy-momentum four vector 

n 

P» = Y,P»s (9) 

s=l 

and with its help define our four vector as 

u, = P,(-P T P T )\. (10) 

We immediately see that when the center of mass of our n-body 
system is at rest, i.e., Pj = 0, % = 1, 2, 3 our four vector w M takes the 
form Eq. (6) in which the wave equation (8) becomes 



r°£po s + £ r °(7 s -P* 



^ = o (n; 



s=l s=l 

where bold face letters mean three dimensional vectors and 



r° = IK, r2 = (7?)- 1 r°- (12) 



r=l 
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Multiplying Eq. (11) by T° and making use of Eqs. (5), (9), (12) 
we obtain 

n 

[-P o + ^(a s -p s + mft)]^ = (13) 

where we used a metric in which P Q = —P° and the latter is the zero 
component of the four vector P M , i.e., the total energy of the system. 

So far we have obtained a Poincare invariant wave equation for a 
system of non-interacting particles which in the frame of reference in 
which P = takes the form (13). 

We wish now to consider interactions and for simplicity we will 
consider them to depend only on the relative coordinates 

x 8 * = x^ s - x,*, /I = 0, 1, 2, 3 (14) 

We want that these relative coordinates become purely spatial ones 
in the frame of reference where the total momentum P = and this is 
easily achieved through the use of the four vector of (6) by defining 

4^<-«vk (is) 

because when takes the form (6) the x s * becomes 

•^o* = ^ ' X T = x, is ~ x it) i = 1? 2, 3 (16) 

As we want our potential to be Poincare invariant it is sufficient to 
make it a function of x^x 5 '*'^ where repeated index /j, are summed over 
their values \i = 0,1,2,3. Thus we restrict ourselves to potentials of 
the form 

V = V(xp s,t '») (17) 

which in the frame of reference where P = becomes V(x st -x st ) with 
the bold face letters indicating spatial relative vectors, i.e. 

x s * = x s -x t (18) 
7 



We are now in a position to write of Poincare invariant wave equa- 
tion for as of n particle of the same mass and of spin | as 

E r * • Wpv + m ) + mrf*"")]* = o (19) 

■8-1 J 

which in the frame of reference in which the center of mass of the system 
is at rest, i.e., P = leads to the Hamiltonian equation 



i) = Eip (20) 



-s=l ^ 

where we replaced P° by the total energy E. So we can denote the 
square bracket in (20) as our Hamiltonian. 

Equation (20) is not the end of our story because even for the one 
particle case it involves both positive and negative energies and it is 
the former ones that will be of interest to us. Fortunately Foldy and 
Wouthuysen [11] gave us a procedure to separate the positive and neg- 
ative energy parts for the one particle case and in the next section we 
proceed to generalize it for the n particle system. 

3 The FW transformation for the many 
body Hamiltonian 

To be able to derive the FW transformed Hamiltonian for the n body 
system we need to review briefly the corresponding analysis for the one 
body case. 

The Hamiltonian is then [12] 

H = + £ + V (21) 

where the odd part O is a • p, the even one £ = (3m and V is the 
potential being a scalar function of the one variable x. 
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We now follow FW [11] and, in particular, the book of Bjorken and 
Drell [12] which states that it is possible to find a unitary operator 



U = exp(iS) (22) 

that allows us to transform H in a series of inverse powers of the mass 
m corresponding to the positive energy part of H. 

The S in (22) is given as 



S = -lL< + 0' + 0") (23) 
2m 



with 



= .p, O' = ,3^1, O" = -^f (24) 



and the new positive energy Hamiltonian H' is now given as 

(p x E) - (E x p) 



+ 4m aS 



(25) 



where one makes the expansion only through order (kinetic energy 2 /m 3 ) 
and [(energy) (field energy)]/m 2 . The field strength E is given by 

E = - W (26) 

and s is the spin of the particle s = —i/Aa. x a. 

Once we have the above well known results we pass to obtaining the 
corresponding ones for the n-body problem 

We start with the two body case where 

H = H 1 + H 2 + V(x 1 ,x 2 ) (27) 

with 

H s = a s ■ p s + (3 s m, s = l,2 (28) 
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Rewriting H as 

H=[H 1 + V(x 1 ,x 2 )]+H 2 (29) 

we can first apply to it the unitary transformation 

E/i = exp(iS'i) (30) 

where Si is given by 5 of (23) where /3, ct, p have the index 1 i.e. Pi, oti, pi 
to indicate its dependence on particle 1. As Si only depends on the 
variables of particle 1 it is clear that it does not affect H 2 in (29) and 
for (Hi + V) it gives the result of (25) with an index 1 for the variable 
i.e. 

exp(iSi)H exp(-iSi) = H x + (H 2 + V) (31) 
where Hi is given by (25) with index one for all the variables. 

We now apply the unitary transformation of the second particle i.e. 

U 2 = exp(iS 2 ) (32) 

where £2 is given by S of (23) where (3, x, p have the index 2 i. e.[3 2) a 2 , p 2 - 
It is clear that the unitary transformation U 2 has no effect on Hi as the 
only terms it could affect are (V\V/&m 2 ) and S • (p x E) — (E x p) 
and this would be of higher order that the ones we accept. We have 
then that U 2 only acts on (H 2 + V) giving us 



U 2 (H 2 + V)U\ = H 2 + V (33) 

where operator H 2 is the one given in (25) where all the variables have 
index 2. Thus if consider U 2 Ui as our unitary transformation we have 

U 2 Ui (Hi + H 2 + V)U\U ] 2 =H 1 + H 2 + V (34) 

The procedure for the two particle problem immediately suggests 
that for the n-body case where V r (xi,x 2 , . . . x„) we can carry out the 
transformation 
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H' = U n U n -! ■ ■ ■ U 2 U 1 HUll4 ■■■Ul = H 1 + H 2 + --- + H n + V (35) 
where 

H t = A (rn + i - + ^s, • (p t xE ( -E t x p t ) 
+ s^V?\/, f = l,2,...n 

We must still keep in mind that in our equation (33) we have to take 
into account that the total momentum P = 0, which we can achieve by 
passing from our coordinate system to the Hamilton- Jacobi one as will 
be indicated in the examples discussed in the following section. 

Thus we extend the approximate FW transformation to the case of 
multi particle relativistic wave equations. For exact FW transforma- 
tions for one- and two-particle systems see refs. [13]. 



4 The two-body problem 



In the case of two particles when the total momentum P = we have 
that 

Pi + P2 = 0, or pi = -p 2 =p (37) 

and we shall denote the corresponding relative coordinate vector as 
r = Yi — r 2 . Besides, for simplicity, we shall take the potential V as of 

2 2 

harmonic oscillator one, i.e., V = mu) A r . 

From (36) and using the notation p, r introduced in the above we 
obtain that 

(p x E) - (E x p) 



4m 2 



Si + S 2 



^ 2 v 



(38) 



as 



Ei = E 2 = E = -VV = 



mw 2 r 



and V\V = V\V = 3 



(39) 
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2 2 

when V = mu) A r . Introducing then the total orbital angular momentum 
L = r x p and using (39) we reduce H' to the form 

H> = (ft + ft) (m + £ - £) + + ■ L + zf- (40) 

V 2m 8m 3 / 4 4m 8m 

where 

S = Sl + s 2 (41) 
is the total spin vector and so the spin values can only be 1 or 0. 

We note also that in Eq. (40) besides terms familiar in non-relativistic 
quantum mechanics there are the matrices fli and /3 2 - They can be cho- 
sen in the form of the direct products [10] 



Pi 





where I 2 is the 2x2 unit matrix, and so 



-h 



(43) 



when I 4 and 0g are the 4x4 unit matrix and 8x8 zero matrix corre- 
spondingly. 

Our only interest in the positive energy part of the wave function 
which implies that (5\ + (3% should be replaces by 2 and thus finally we 
have to deal with the expression 

J 2 \ fp 2 mu 2 r 2 J 2 „ T \ p 4 
2m + 3— + - + — ;— + -;— S-L -/-T 
8m/ \m 4 4m / Am 6 



H' 



(44) 



As the second parenthesis correspond to an harmonic oscillator with 
spin-orbit coupling whose eigenfunctions and eigenvalues are very well 
known we can use the former as a complete basis to convert H' into a 
numerical matrix. 
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Before proceeding to give a procedure to calculate the eigenvalues 
for H' of (44) it is convenient to make the canonical transformation 



77 = \ — P> p = \l^ r ( 45 ) 



to write H' as 



where in the second parenthesis we have the oscillator hamiltonian of 
unit frequency and the spin-orbit term remains uncharged as r x p = 
p x 7r while p 4 is replaced by |m 2 u; 2 7r 4 . 

To convert H' into a numerical matrix we can use the states of the 
harmonic oscillator and states of squared total angular momentum J 2 , 
total orbital momentum L 2 and total spin S* 2 , i.e., 

\nl, (^^jS; j,m >=Y, < lf^, S(r\jm > \nlfi > |Q-^cr> (47) 

where J 2 , L 2 and S* 2 commute with the Hamiltonian (46) and so are 
integrals of motion. The kets \nlfj, > are those of the harmonic oscillator 
of unit frequency [14] and | (jj^JSa > are those of the total spin. 

The numerical matrix we want to determine has then the elements 

< n% ^S;j,m\H'\nl, (±£)S,j,m> 
= (2m+ 3 £ + uj(2n + l + l)+ + 1) - 1(1 + 1) - s(s + 1)])<W 

< n'lWnl > 

(48) 

The last term in (48) can then be defined starting with the relation 
< n'/|7r 4 |n/ >=< n'l\iT 2 \n"l >< n"l\ir 2 \n'l > and using the expression 
for < n'l\ir 2 \nl > given in p. 7, Eq. (3.11) of reference [14]. Thus we 
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get that 

< n'l\ii 4 \nl >= yjn(n - l)(n + l + \)(n + l - \)5 n > „_ 2 
+ (4n + 21 + l)^n{n + l + \)5 n , n _i 
+[(2n + l + l){2n + l + \) + 2n(n + I + \)\5 n . n (49) 

+ {An + 21 + 5)^{n + l){n + l + \)5 n , n+1 
+ ^(n + l)(n + 2)(n + l + \)(n + l + |)5 n , n+2 

We note that as we are using units in which h = c = 1, the dimen- 
sionless term u/m becomes in the c.g.s. units 

uj huj 



m mc 2 



(50) 



and can be treated as a small parameter. In following section we com- 
pare the results of our analysis with the experimental spectrum of the 
bottonium and charmonium masses. 



5 Energy spectrum of the two body prob- 
lem 

We start with the eigenvalue problem 

#V = E nl ^ (51) 

for the two-body hamiltonian (48). 

To find the related energy spectrum we are supposed to diagonalize 
matrix (49) which can be done using numerical methods. Moreover, 
for sufficiently small coupling constant (50) it is possible to apply the 
standard perturbation theory and express the eigenvalues of H' (48) in 
power series of v = u/m: 

o 2 

E n i-2m- — = E° nl + E l nl + E 2 nl + ... (52) 
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where the non-perturbed levels E® t are linear in v, 

E° nl = mu(2n + l + ^j (53) 

while E x nl and E 2 t are quadratic and cubic in v respectively. Moreover, 

the first and second perturbations of energy spectrum can be expressed 

via the elements of the perturbing matrix < n'V j'\K\jnl >= f^[j(j + 

1) - 1(1 + 1) - s{s + l)]6 nn ,6 a > ~mi< nT\n*\rd > 5 jf as follows [1] 

^1 , , -rri , ^9 <n'l'j\K\n'l'j> 2 , 
E 1 nlj =<nlj\K\nlj>, E 2 nlj = £ J' V ( 54 ) 

Then, using (49), (52) and (53) we obtain 
E nl = 2m+?£+u;(2n + l + l) + £-[j(j + l)-l(l + l)-s(s + l)} 



Am 



16m V 2 / 8m I- 

[(2ra + / + |)(2n + / + |) + 2n(n + Z + §)] 



\-" + ^ + 2 



9n(n + Z + |) + 2^2n + Z + f^2n + Z + 



(55) 

Formula (55) describes the spectrum of relativistic two body systems 
with the harmonic oscillator potential. We compare it with the energy 
spectrum of two quark systems (mesons) and find that it presents a 
rather realistic qualitative distribution of bottonium and charmonium 
masses. For example, setting m = 4.7 GEV (i.e., supposing m be equal 
to the bottonium quark mass) and choosing the dimensionless coupling 
constant v — 0, 19 we obtain from (55) the following values for the mass 
differences of the bottonium states: £f,o — T(IS') = 0.411 GEV, ^i(lP) — 
T(1S) = 0.423 GEV,£ b2 {lP) - T(1S) = 0.443 GEV while the experi- 
mental data are [15] 0,400, 0,432 and 0,453 GEV respectively. 

To obtain a better agreement with experimental data we consider 
anharmonic oscillator potential 

V = ^ + V>, V' = - am ^ rA (56) 

where a is a dimensionless interaction constant. In addition, we take 
into account relativistic corrections up to order u 3 for the approxi- 
mate Hamiltonian H' which needs continuation of the FW reduction. 
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We omit the routine calculations which are analogous to ones given 
in Section 3 and present the resulting transformed Hamiltonian using 
variables (45): 

^ = 2- + f^ + |(- 2 + P 2 )+£S.L-^(^ + a/) 

+ 32^ (t + ( 2 " ha )P 2 ~ ^ + 2 / ?2 ) S • L - (S • tt) 2 + (S 2 - 8)tt 2 ) 

(57) 

We need only diagonal matrix elements for terms of order u 3 /m 2 , 
placed at the second line of equation (57). They can be easily found 
starting with matrix representation for p 2 and 7T 2 given in page 7, Eq. 
(3.11) of reference [14] and the representation for S • 7r in the spherical 
spinor basis given in pages 422-423 of reference [16] : 



< nl\p 2 \nl >=< nl\n 2 \nl >= 2n + I 



3 
2> 



< nl\p 4 \nl >= (2n + l + |)(2n + / + §)+ 2n(n + I + 
< nl\7i 6 \nl >= 4{n + l + \){2n + 1 + 1) 
+ {2n + 1 + 1) (2(n + l)(n + I + |) + (2n + I + |) 2 ) , 
< nlj\(S ■ n) 2 \nlj >= (l - \{j - I) 2 + ^y) < nl\n 2 \nl 

and the related energy values E n ij are described by the following for- 
mula 



> 



E nlj = 2m+ + u{2n + l + \) + ^ + 1) - /(/ + 1)] 

_(i+aV T( 2 n + 1 + 1) (2n + l + l)+2n (n + l + \) 



16m 

,.,3 



256m' 



7 (2n + I + |) (1 - a) 2 (2 (2n + l + l) 2 + 4n {2n + l + \ 



+ (1 + a) 2 Un + 1) (n + I + f ) - \ [2n + I + \ 



32m 2 

a; 3 
64m 2 



(2n + / + |) (j(j + l)-l(l + l) + 5-S 2 ) 
(4n(n + / + i)(2n + / + l) + (2n + / + |) (2(n + l)(n + / + |) 
+ (2n + / + |) 2 )) 

+ ^,(2n + / + |) [3Z(Z + 1) - 3j(j + 1) + (J - 2 + £b + 5S 2 - 14 

(58) 

We compare the spectrum (58) with the experimental bottonium 
and charmonium mass spectra [15]. We set m = A.7GEV for the bot- 
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tonium and m = 1.4 for charmonium cases and use MAPLE software to 
find the coupling constants uj and a which correspond to the minimal 
deviation of the spectrum (58) from experimental data. The results of 
our investigations are given in the following tables. 

Table 1. Experimental spectra of the bottonium and our model 
results (in GEV), m = 4.7, u: = 0.378, a = 4.7 

Experimental Theoretical 



1 = 0,3 = 


1 






T(2S) - 


-T{1S) 


0.563 


0.574 


T{3S) - 


-T{1S) 


0.895 


0.963 


T{4S) - 


-T(1S) 


1.119 


1.142 


1 = 1,3 = 









6o(iP) 


-T(LS) 


0.400 


0.310 


6o(2P) 


-T(IS') 


0.772 


0.797 


1 = 1,3 = 


1 






61(1^) 


-T(IS') 


0.432 


0.316 


6i(2P) 


- T(1S) 


0.794 


0.801 


1 = 1,3 = 


2 








-T(IS') 


0.453 


0.326 


61 (2P) 


-T(IS') 


0.808 


0.808 
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Table 2. Experimental spectra of the charmonium and our model 
results (in GEV), u = 0.45, a = 1.2, m = 1.4 



Experimental Theoretical 



l = 0,j = l,s = 



r)' c (2S) - r) c (lS) 0.603 

l = 0,j = 1,8 = 1 



0.604 



V(2S) - J/V(1S) 0.589 
V(3S) - J/tf (IS) 0.943 
/ = l,j = 0,s = l 



0.589 



0.925 



U(1P) - J/*(1S) 0.308 
/ = l,j = l,s = l 



0.372 



&i(2P) - J/*(1S) 0.413 
/ = l,j = 2,s = l 



0.390 



e c2 (3P) - J/V(1S) 0.459 



0.422 



We see that in average the mass spectrum predicted by our very sim- 
ple model is in rather good accordance with the experimental data. It is 
possible to obtain a better agreement with experimental data changing 
the quark masses by effective ones which are additional free parameters. 

6 Conclusion 

We propose a relativistic Hamiltonian for an n-body problem in which 
all the masses are equal and all spins take value 1/2. Discussing the 
problem in the frame of reference in which the total momentum P = 0, 
we were able to extend the FW transformation to n-body case and 
separate the positive energy part of the Hamiltonian. Examples of two 
body systems are discussed in more detail. 

The proposed approach admits a straight forward generalization to 
the case of particles with different masses and spins and is valid for 
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more general form of the interaction potential. In the present paper 
we discuss only the main ideas and demonstrate it effectiveness using 
the simplest interaction model. Nevertheless, even this very straight 
forward model predicts a rather realistic bottonium and charmonium 
spectra presented in the Appendix. It looks rather curiously that start- 
ing with a an interaction potential which is not well grounded physically 
and using only two free parameters uj/m and a it is possible to obtain 
a good qualitative and also relatively good quantitative description of 
bottonium and charmonium masses. 
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